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ABSTRACT 

A general treatment of disk star formation is developed from a dissipative multi-phase model, 
with the dominant dissipation due to cloud collisions. The Schmidt-Kennicutt law emerges naturally 
for star- forming disks and starbursts. We predict that there should be an inverse correlation between 
TuUy-Fisher law and Schmidt-Kennicutt law residuals. The model is extended to include a multi-phase 
treatment of supernova feedback that leads to a turbulent pressure-regulated generalization of the 
star formation law and is applicable to gas-rich starbursts. Enhanced pressure, as expected in merger- 
induced star formation, enhances star formation efficiency. An upper limit is derived for the disk star 
formation rate in starbursts that depends on the ratio of global ISM to cloud pressures. We extend 
these considerations to the case where the interstellar gas pressure in the inner galaxy is dominated 
by outflows from a central AGN. During massive spheroid formation, AGN-driven winds trigger star 
formation, resulting in enhanced supernova feedback and outflows. The outflows are comparable to 
the AGN-boosted star formation rate and saturate in the super-Eddington limit. Downsizing of both 
SMBH and spheroids is a consequence of AGN-driven positive feedback. Bondi accretion feeds the 
central black hole with a speciflc accretion rate that is proportional to the black hole mass. AGN- 
enhanccd star formation is mediated by turbulent pressure and relates spheroid star formation rate 
to black hole accretion rate. The relation between black hole mass and spheroid velocity dispersion 
has a coefficient (Salpeter time to gas consumption time ratio) that provides an arrow of time. Highly 
efficient, AGN-boosted star formation can occur at high redshift. 

Subject headings: galaxies-disk: galaxies: elliptical - galaxies: active galactic nuclei-galaxies: evolu- 
tion -stars: formation 



(N 
> 
O 
00 

o 

o\ 
o 



1. INTRODUCTION 

Supernova feedback is considered to be a crucial el- 
ement for negative feedback in star formation in disk 
galaxies. The star formation history in massive spheroids 
requires, according to the prevalent view, negative feed- 
back from AGN. Whether this is sufficient to explain the 
observed downsizing is far from clear. Here we reassess 
the Schmidt-Kennicutt (SK) star formation law and de- 
velop a simple multi-phase model i n terms of the porosity 
formalism applied to disk galaxies (|Silkll200i[ l.We extend 
the model to incorporate AGN-triggered star formation 
and provide an application to spheroid formation and 
ultraluminous starbursts. 

A cloud collision mod e l of th e SK law has been previ- 
ously presented by iTanI ()1999[ ). who uses galactic shear 
to compute the cloud collision rate. One advantage of 
this approach is that it provides a natural explanation 
for the low star formation rates observed in the outer 
parts of disk galaxies and complements an alternative 
explanation which appeals to UV background radiation- 
control l ed H2 suppression in the dust-deprived outer disk 
ISchavd (|2004j) . We provide a simplified reformulation be- 
low, that we will apply in the context of a multi-phase 
medium to incorporate star formation and supernova 
feedback (Section 2). In Section 3, we explore regula- 
tion of star formation by turbulent pressure and set an 
upper limit on the disk surface brightness due to star 
formation. Section 4 builds on the AGN feedback model 
[Silk (2005) and applies AGN triggering to star formation 
in protospheroids. Scaling laws are derived for the black 
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hole growth rate and the star formation rate. Downsiz- 
ing of both super-massive black holes and stellar mass is 
found to be a natural consequence of Bondi accretion-fed 
black holes and AGN-induced star formation 



DISK STAR FORMATION RATE: 
MODEL 



CLOUD COLLISION 



Consider cloud collisions in the disk as a trigger of star 
formation. Cloud formation and collisions are driven by 
the non-axisymmetric gravitational instability of a cold 
self-gravitating gas-rich disk. Let a typical cloud have 
pressure pd and surface density Sc/- We expect star- 
forming clouds to be marginally self-gravitating and also 
to be confined by ambient gas pressure. Clouds form this 
way, and may be maintained if the cloud covering factor 
is of order unity, this condition guaranteeing that colli- 
sions occur on a local dynamical time-scale. If the clouds 
are strongly bound, it is difhcult to avoid a short lifetime, 
collapse and star formation. Our description is a statis- 
tical one where we are assuming a steady state ensemble 
of clouds although the clouds are being formed and re- 
formed all the time in competition with cloud destruction 
and dispersal processes such as star formation and col- 
lisions. For typical parameters in the Galaxy, the disk 
crossing time normal to the disk and the cloud lifetimes 
are similar both of order lOMyr, although we assume in 
general a statistically steady state cloud population in 
this analysis. 

The following relations then apply: 

Pg = Pg'^l = TrCEgEtot, (1) 
assuming equal scale-heights for clouds and stellar mass. 
If the clouds are self-gravitating, then p'^J — ttxGSc/^, 
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where x ^ 10 is an estimate of the pressure en- 
hancement due to self-gravity of interstellar clouds. 
We redefine = XPci, and can now write Sc; = 
{Vci/PgY^'^iJ^tot^gY^'^ ■ The covering factor S^-i of clouds 
in the disk is directly inferred to be Sd = {^g/'Sci)fch 
where fd is the gas fraction in clouds. We rewrite this 

as Sd = fdiPg/Pci)^^'^{,^g/^tot)^^^ ■ Here Eg is the total 
(cloud plus diffuse) gas surface mass density. The cloud 
collision time-scale is tcoii = {TidH) / {Yigfd<Jg), where 
the scale height H^^ = (7rGStot)/crg, and dg is the cloud 
velocity dispersion. The collision time can also be ex- 
pressed as tcoii = Sd^^tcross with tcross = H/ag, which 

becomes tcoii = f^i^{pd/PgY''^{'£'tot/^gf^'^{H/ag). 
More gener ally, inclusion o f more realistic 3D cloud kine- 
matics (cf. iTasker fc TanI (|2008) ) yields correction fac- 
tors of order unity. 

We now assume the disk star formation rate is 
self-regulated by supernova feedback which drives the 
cloud velocity dispersio n. While this ass umption 
has a long history (cf. iFirmani fc Tutukov I ([1992)), 
it remains controversial. Numerical simulations cer- 
tainly demonstrate that supernovae provide negative 
feedba ck into star-forming clouds by driving t urbu- 
lence (Joung fc Mac Low* '2OO6'; 'Tasker & Brvan' '20061; 
Kovama fc Ostrike r 2008a,b; Kim fc Ostriker 2007; 
Joung et"aLl i2008f) . Turbulent pressure plays an im- 
portant role in regulating star formation, via control- 
ling theporosity of supernova remnant-driven bubbles 
(Silk 2001) as well as the molecular hydrogen fraction 
(Blitz fc Rosolowskvl |2006[) . At the same time, global 
shear also plays a role in controlling cloud peculiar 
veloc ities, especially for massive clouds (jGammie et al.l 
|1991|) . Since global gravitational instabilities ultimately 
drive cloud formation, and hence control star formation, 
the common origin of competitive turbulence drivers 
means that effects of shear and supernovae in self- 
regulating cloud turb ulence are not easily sepa rated in 2- 
dimensional models (jShettv fc Ostrikeill2008f ). However 
fully three-dimensional high resolution models of self- 
consistent star-forming disks embedded within dark ha- 
los demonstrate that non-axisymmctric gravitational in- 
stabilities dominate the observed turbulence of ~ lOkm/s 
at low star formation rates, but that supernova feed- 
back will be important via the intermediary of the hot 
gaseous phase at a star fo rmation rate in excess of 
10::iMQkpc~2yr-i (jAgertz et al.l 120091 : ITamburro et all 
[2OO9I) . 

Let TOsAT be the mass in stars formed in order to result 
in a Type II supernova. This is just a function of the 
adopted IMF. Momentum balance gives 

T,*{EsN/{msNVc)) = fc^gO-g/tcoU- (2) 

Here fc is the cloud volume filling factor, which can be 
expressed in terms of porosity Q as fc = . Also, Esn 
is the kinetic energy of a SNe II and Vc is the velocity 
at the onset of strong cooling of the SNe II remnant. 
Canonical numbers used throughout are msN = 15OM0 
(for a Chabrier IMF) and Vc = 400kms~^. 
We can rewrite the star formation rate per unit volume 

as 

P* = eSNfcfcl\/GPgPg (3) 

with esN = {msNVc<Jg)EsM~^{,Pg/pdY^'^- This formu- 
lation is commonly used as a star formation rate pre- 



scription in semi-analytical modeling of galaxy forma- 
tion. It may be more relevant to rewrite this formulation 
for disks: 

= fcfdG{nll,o,f' h]^NVc)\ [^^/^^^ji 

= esNfcfcl^/Tg{R/H)^^gas^ (4) 

where T,tot — Sg -I- S*. Here the disk gas fraction is 
fg ~ 0.1, and we use the disk scale- height -to-radius 
relation H/R = [ag/vrY ^'^^ a disk rotating at Vr with 
Sl^ = GYjtot /R- Remarkably, although the preceding for- 
mula ignores the multi-phase nature of the interstellar 
medium and the possibility of gas outfiows (see below), 
one nevertheless manages to fit the Schmidt-Kennicutt 
relation. 

Wc write the observed Schmidt-Kennicutt (SK) law as 

3 /2 

S* = CsK^g 1 and obtain the SK law coefficient 

CsK = TT'^^Gfcfd{msNVc/EsN)[Pg/Pcl]--^]i!. (5) 

Inserting typical parameter values, we find that 

This demonstrates that we get the correct normalization 
at, say, 3 kpc, the scale length of the molecular gas in 
the Milky Way, where the scale-height is around 100 pc, 
the gas fraction is around 0.2, and the molecular gas cov- 
ering fraction around 30%. The observed star formation 
efficiency in inner spiral disks is found to be fairly ro- 
bust and for H2 a lone amounts to 5.25 ± 2.510~^°yr~^ 
(|Lerov et al.ll2008D . 

This compares well with the Kennicutt law, both lo- 
cally and at z ~ 2 in shap e (!]» ^ Spas) and in normal- 
ization (for esN « 0.02) (jBouche et al.il2007[ ). For the 
luminous star bursts at z 2, the turbulence is enhanced 
(fjg ~ 40km/s), but the scale height is thickened. One 
reason is that esN oc ag and {R/HY^'^ cx l/dg for disks 
with varying amounts of turbulence, as might be induced 
by minor mergers. If the covering factor increases, it is 
not obvious if the star formation rate in a cloud colli- 
sion model would increase. To lowest order, these effects 
all cancel at fixed Tttot-, and we can hence understand 
how starbursts remain on the local Schmidt-Kennicutt 
law. Supernova feedback effectively keeps star formation 
inefficient. 

Of course, we need to better understand how starbursts 
satisfy the same scaling law as quiescent disks. One hint 
is that while the gas velocity dispersion may vary in star- 
bursts depending on the merging history, Sjot satisfies 
Freeman's law and is approximately constant for star- 
forming disk galaxies. The observed dispersion in the 
Schmidt-Kennicutt law may arise from the dispersion in 
total surface density and molecular as well as total gas 
fraction fg. 

2.1. TuUy-Fisher relation 

The TuUy-Fisher relation is also controlled by disk sur- 
face density. We use the empirical I-band TuUy-Fisher 
(TF) relation: = Ctf v'^, with a « 4 in the K- 
band ({Masters et al.l |2008[ ) and Vr the maximum rota- 
tion velocity, and where the virial theorem requires that 
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Ctf = (3/47r)G-2l]-^i(L,/Aftot)- We find using equa- 
tion ^ that 

(7) 

We infer that the Schmidt-Kennicutt law residuals 
should anti-correlate with the TuUy-Fisher law residuals. 
The TuUy-Fisher normalization is correct, by assump- 
tion: what is new is the predicted inverse correlation 
between SK and TF law residuals. 

2.2. Gas- dominated disks 

The global star formation law can be applied to re- 
gions that are gas-dominated. When gas dominates 
the self-gravity, the cloud collision model suggests that 

S* cx Sg'^/^Stot"^^^ ~ Sgos^, and the KS law steepens. 
There are indications of such a steepening in several en- 
vironments. 

1. The extended H I spiral structure in NGC 6946 
(jBoomsmal l2007| ) shows that global gravitational 
instability is not a sufficient condition for form- 
ing stars. In the case of M83, the HI disk ex- 
tends to more than twice the optical scale. Deep 
UV imaging reveals very low level star formation 
in the outer HI disk, well below the SK thresh- 
old. The cloud collision model provides a possible 
explanation of these phenomena, although the ob- 
served radial dependence of star formation rate is 
too steep to be ex plained by the simplest models 
()Lerov et al.ll2008f l. 

2. Individual young star complexes in M51 fall on 
the SK law, although with increased dispersion 
(jKennicutt et al. 2007) and a slightly steeper slope. 

3. Damped Lyman a lpha systems at z ~ 2 
() Wolfe fc Chenl[2006l ) underproduce stars by up to 
a factor of 10 in star formation rate as predicted 
from the SK law. 

4. Steepening also occurs in the inner regions of disks 
at extreme star formation rates. This is found 
i n intensely star -fo rming galaxies at high redshif^t 
()Gao et al.'ll2007l ). iKrumholz fc Thompson! ()2007l ) 
account for the linear relation found for the local 
HCN data in terms of the critical density for exci- 
tation of the HCN transition in dense gas, that ef- 
fectively samples only the densest molecular clouds 
and thereby bypasses the sensitivity to dynamical 
time-scale. 

Steepening in a cloud collision model is not a unique ex- 
planation for any of these phenomena. For exam ple, the 
outer parts of disks are more thermally stable (jSchavel 
|2004[ ). and the star formation rate in DLAs could be 
suppressed because of the low H2 content due to a com- 
bination of a low dust content plus a high radiation field. 

3. PRESSURE-REGULATED STAR FORMATION AND 
STARBURSTS 

Turbulent pressure-regulated star formation is espe- 
cially likely to be important in starbursts. In disks, 
atomic cooling provides an effective thermostat for the 



turbulent velocity dispersion. Feedback operates via 
the hot phase venting into the halo. Gas may cool 
and fall back into the disk, as in the galactic foun- 
tain model, or escape in a wind, as happens for dwarf 
starburst galaxies. The simple porosity des cription o f 
supernova feedback in a multi-phase ISM l200lh 
provides an expression for the star formation rate in 
which porosity-driven turbulence is the controlling fac- 
tor: jO, = QmsNi'^'^ l'iR^ta)~'^ , where the shell reaches 
a final size Ra before break-up, determined by the ambi- 
ent pressure at expansi on time tg. The sh ell evolution is 
generally described by (jCioffi et al.lll988[ ) 

t^t,El{'\-^l\vJvY^I' (8) 

and 

R = R,El{'n-^l\tmf''', (9) 

where vq = 413km/s, Rq = 14pc, to — 1.3 x 
lO^yr. Here cooling becomes significant at shell veloc- 
ity Vc = 413i?5(^ng^'*A'^/*kms~^ where the cooling time- 
scale within a SN-driven shell moving at velocity Vc is 

tc = Vc/Xp, 

A-i = 3m3/2fci/2j.i/2/A^^^ (10) 

and Aeff{T) is the effective cooling rate ( ^ f^^^ over 
the relevant temperature range 10"* < T < lO^K associ- 
ated with cooling shock velocities < lOOfcms"^ ). 

The SNR expansion is limited by the ambient turbulent 
pressure to be PgCg, and we identify Va (shell velocity at 
time ta) with (Tg. We obtain 

- Q^gPg{ag/af,dy^/\ (11) 

where 

a,,:. = (coGV2„^3/2„19/7^62/49^-^)V19„-1/14 (13) 

« 20n-'^'^mf^%,E",fkn,/s (13) 

and Co = ^i?oio- 

The dependence of star formation rate on turbulent ve- 
locity is reminiscent of Barnes's model for star formatio n 
in the Mice galaxies, an ongoing merger (|Barnesll2004f) . 
A turbulence prescription is required to reproduce the 
observed spatially extended stellar distribution, which is 
inconsistent with the simple density-dependent Schmidt- 
Kennicutt law. 

We apply the cloud collision model of Section 2 to com- 
pute Q. The star formation rates derived via porosity 
and via cloud collisions can be set equal. Comparison 
with the star formation rate derived from cloud collisions 
yields 

Q=(^-^Y%M^^f^(^Y\ (14) 

then using fc = 6^*3, we find 

Qe"^ = fclesNJ^, l^^y^' . (15) 

This is an explicit expression for the porosity as a func- 
tion of the turbulent velocity. The star formation effi- 
ciency is evaluated here at the fiducial velocity dispersion 
so that: 

f-SN.Jid = {msNVcCrfid)EsN^^{Pg/Pcl)^^'^ ■ (16) 
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There are two regimes: Q << 1 and Q >> 1, where 
approximate solutions can be found and the real solution 
joins them smoothly. For canonical values, with pd ~ Pg 
and 17 g ~ (Jfid the right hand side of the above equation 
for Q is small (due to the 1% efficiency of star formation 
calculated previously) and is of order esN.jid ~ 10~^ and 
thus 

Q ^ fclesN.J^d (^^) . (17) 

For completeness we give here the case where Q >> 1, 
namely 



Q-ln 



fcl^SNJid 



12/7' 



(18) 



An approximate formula encompassing both regimes 
Q » 1 and Q « 1 is 



Q-ln 



1 + fcl^SNJid 



(Jfid 



12/r 



(19) 



Generally, we find that the star formation efficiency is 

- n-7/12„-7Q/12/ n19/12„7/12 ,„„^ 

A desirable feature is that the star formation rate van- 
ishes at very large Q. and becomes larger for small Q. 
There is no minimum in the star formation efhciency- 
as a function of Q but the above features suggest that 
Q asymptotically becomes constant, of order unity, and 
self-regulation occurs since as Q exceeds unity it depends 
only logarithmically on the velocity dispersion. One 
can better understand why starbursts lie on the same 
Schmidt-Kennicutt law if we assume that local physics 
specifies the gas fraction converted into s t ars, i n effect 
esN, as in the model of iKrumholz et all (|2006D . This 
is plausible for individual molecular cloud complexes. A 
constant star formation fraction (equivalently, efficiency) 
is also expected globally in quiescent disks. Since Vturb 
self-regulates at ag ^ lOkm/s, then if porosity also self- 
regulates, the efhciency or fraction of gas converted into 
stars per dynamical time is constant and small. Then, 
the higher turbulence in a starburst means that the 
porosity is low. In turn, low porosity guarantees inef- 
ficient feedback and runaway star formation. 

In merger-driven starbursts, the porosity is small since 
esN ~ CTg, whence Q (x ag Small Q suggests a nu- 
clear starburst, whereas large Q regulates global feedback 
and the star formation rate in a disk. This is complicated 
by the dependence of eg n on Vturb which is compensated 
by the increase in scale-height with enhanced turbulence. 

In a quiescent star-forming disk galaxy, we might ex- 
pect the porosity to self-regulate and be of order unity. 
This is the case, for example for the Milky Way, where 
the supernova feedback regulates the global star for- 
mation rate. Even in this case, the star formation is 
not monotonically decreasing with time, as the Schmidt- 
Kennicutt l aw w ould suggest. Numerical simulations 
(|Slvz et al.l 120051 ) suggest that mini-starbursts occur 
stochastically, on a scale of order 1 kpc, with the mean 
global value decreasing as the gas supply is reduced. Ob- 
servational evidence for a non-monotonic star formation 
history in the solar neighborh ood comes from surveys o f 
chromospheric age indicators (|Rocha-Pinto et ani2000f ). 



Evidence for a series of starbursts is found in the inner 
disk star- forming regions of spiral galaxies (jAllard et al.l 
[20061) . 

Another aspect is the extreme pressure induced by gas 
dissipation. This must play a role in regulating star for- 
ma tion. An explicit cas e for p ressure regulation is made 
by iBlitz fc Rosolowskvl (|2006f ) who propose a modified 
Kennicutt law: oc Egp°-^. Our SN-regulated law us- 
ing equation ^ is 

= G^'^{msNVc/EsN)[PglPclY'^Pg^'^^g, (21) 

becoming 

= TT^^/'^msMVc{EsN^ciy^Pg'i^g- (22) 

This is close to a pressure-regulated star formation law, if 
all star-forming clouds have a threshold column density. 

3.1. Outflows 

If the porosity is high, as may happen transiently in 
starbursts, we assume that disk outflows occur. These 
may be winds from dwarf galaxies or fountains in the 
case of more massive disks. Numerical simulations of 
star formation in the multi-phase interstellar medium of 
a disk galaxy are able to mo del the disk outflows and 
global star formation history (jTasker fc Brvanll2006f ). It 
is useful however to provide an analytic formulation. Let 
Jl be the hot gas loading factor. It is measured to b e 
around 10 for the outflow from NGC 1569 (|Martinll2005l) . 

Now the outflow from the disk is 



Mont = (1 - e-^)fLM* ~ Q/lM* 



(23) 



This tells us that Mout/M^ = IlQ ^ ^^g^'' ■ Low porosity 
suppresses outflows so that outflows are suppressed in 
massive potential wells. 
We can also express the outflow rate as 

Mout ~ Q^fL{aglaf,df-'fl'^MglU. (24) 

This shows that outflows could indeed occur from mas- 
sive potential wells if porosity can somehow be main- 
tained. We argue below that AGN-triggered star forma- 
tion fulfills this role. 

Outfiows are important for dwarf galaxies, but are seen 
to be quenched in deep potential wells as well as at ex- 
treme porosity. If Ug is high, the porosity is low, feedback 
is suppressed and supernova-driven winds are quenched. 
Even if Ug is low, outfiows may be suppressed if the gas 
is dense, leading to low porosity. Dwarf galaxy outfiows 
play an important role in IGM enrichment at high red- 
shift. The mass ejected during dwarf formation is com- 
parable to the mass retained in stars formed. Nearby 
starbursts display this trend, suggesting that the effect 
may be generic to powerful starbursts in dwarf galaxies. 

3.2. Upper limit on the disk star formation rate 

iMeurer et all (|1997[ ) and iHathi et al.l (|2008f l report a 
bolometric upper limit on the disk surface brightness in 
starbursts of 2 x lO^Lgkpc'^ over 0.1 ^ Re ^ lOkpc. 
We interpret this upper limit on disk luminosity for our 
model galaxy in terms of radiation pressure and mechan- 
ical pressure from limi ting the gas surface density (cf. 
[Thompson et al.l (|2005f ). 
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To avoid lift-ofF via application of the Eddington con- 
dition requires: 
(1) for radiation pressure: 

Si < {TT/2)Gc^tot^g (25) 



S* < cnGfi ^T,tot^g/2, where /; 



10- 



10 c for massive stars. Let us evaluate this expres- 
sion after using the earlier expression for to eliminate 
Eg. We find using equation (|25p that star formation is 
reduced at low gas surface density and lift-off can be 
avoided if 
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(2) For mechanical energy input: we have instead no lift- 
off provided that 
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(27) 



Here the symbol e,„ denotes the efficiency of conversion of 
rest mass energy to mechanical energy in massive stars. 
For the lift-off limits to be reached the left-hand side 
of the above two equations must be greater than unity 
since T,g < Stot. For the radiation case, unity is not nor- 
mally reached on the left hand side unless some of the 
parameters such as ci are signific antly different from their 
normal values. This is what led [Thompson et al.l (120051 ) 
to argue for an AGN luminosity driving this Eddington 
limit for disks. For mechanical winds, the left side is 
boosted over the equivalent radiation driving term by 
a factor of {c/v^Y making it reasonable that such me- 
chanically driven feedback occurs. For completeness and 
because of our poor knowledge of star formation at high 
redshift we will keep both the radiation and mechanical 
estimates in our calculations. This upper limit on Eg 
corresponds to an upper limit on global star formation 
of 45M0kpc~^yr^^. By writing the disk surface bright- 
ness as Ei = ejc^E*, we find using equation (|26p that 
the disk surface brightness satisfies 



Ei < (|)2GcEL 



EsN 



Pc 



(28) 



For mechanical input we have a similar equation using 
equation (P7)l 



Ei < (|)^GcEL 



Es 



N 



msNC^ £71 



Pd 
Pg 



^— . (29) 



We see that the absolute maximum global star formation- 
driven surface brightness, ^L,max, is 



'^L,max = (7r/2)GcE 



2 

tot- 



(30) 



This expression can be derived directly from the pre- 



vious equation by setting Eg 



Note in this 



case all the photon energy is used up to support the 
disk and the photon bubbles can occur with the sig- 
nature of light and dark patches on the disk. This 
suggests saturation of the star formation rate occurs 
at high surface density. Inserting typical numbers, 
Ej < lO^Le kpc-^EiiV-^ooe-iaEt^^t iooQ. Note that 
iKomugi et al.l ()2005l ) report an offset roughly equiva- 
lent to an effective steepen ing of the S K law for ultra- 
luminous starbursts, as do iGao et al.l ()2007| ). We now 



find that for an extreme starburst at the Eddington lu- 
minosity with Thomson scattering, ax, as the dominant 
opacity the relation between surface brightness and mass 
surface density is: 



'E 



(31) 



where rUp is the proton mass. From ((28|) and ((29|) we 
deduce that the total surface density, which is essentially 
all gas in initial situations, satisfies 



> E, = — 



EsN 



(^) 



Pg_ 

Pel 



nip 
(7T 



~ lO^Mg pC"^. 

and similarily for mechanical input: 



(32) 



> E, 



EsN 



Pg_ 

Pel 



lILp 



(33) 

A direct and simple way to derive these limits is to take 
equations ((3T|) and (|26p and find that 



> 



iii,p 



4 X lO^Afopc" 



(34) 



which is a derivation of Fish's Law. The densest galactic 
molecular clouds, trac ed by H2O ma sers, have a simi- 
lar surface density Plu me et al.l ()1997f ). By z ~ 2, there 
is a substantial population of galaxies with star forma- 
tion rates of 5OO-lOOOM0/yr, ULIRGS an d SMGs. A re- 
cent st udy of SMGs at 0.5 arcsec resolution lTacconi et al.l 
(|2008f ) finds that they are gas-rich (molecular gas frac- 
tion ^ 30%) with the gas in compact disks at a density 
of ^ IO^Mq pc~^, and are undergoing major mergers. 

4. APPLICATION TO ACTIVE GALACTIC NUCLEI 

AGN outflows over-pressure the interstellar medium. 
They can deplete the gaseous environment by driving 
a wind. AGN outflows are the principal element in 
semi-analytical modeling of massive ellipticals that helps 
quench recent star formation. The energetics are as fol- 
lows: the specific energy per baryon from supernovae 
is Esn/itT'SN ~ 10~^c^ergs/gm, whereas AGN outflows 
provide ~ 10~^c^ergs/gm per unit spheroid mass, for 
an assumed efficiency of O.lc^ and a SMBH-to-sphcroid 
mass ratio of lO"*^. We argue below that AGN outflows 
have global impact by driving overpressurised cocoons 
into the inhomogeneous ISM. 

Useful formulae are: 
(1) the Eddington luminosity 



LEdd = iTTcGMBHrnp/cTT, 

(2) the Salpeter time-scale 

ts = rjc^MBH / Lsdd = rjaTciATiGmp)'^ , 

and 

(3) the self-regulating feedback mass 



MsRF = fc 



nip ttG^ 



(35) 
(36) 

(37) 



Blow-out by radiation pressure occurs (assuming a ho- 
mogeneous interstellar medium) if L = Lsdd at Mbh = 
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MsRF (|Silk fc Reeslll998l: iFabianlfToMl . This will lead 
to a wind, deplete the baryon reservoir and quench black 
hole growth by gas accretion. With fg ~ 0.1, as ex- 
pected initially in the protogalactic core, this simple re- 
lation fits the mean of the observed relation over 3 or- 
ders of magnitude in black hole mass. Numerical simu- 
lations generally confirm these simplistic estimates. We 
further define a critical AGN luminosity Lsrf for star 
formation-boosted AGN outflow by the Eddington lumi- 
nosity associated with the critical black hole mass Msrf 
that corresponds to the balance between Eddington lu- 
minosity and proto-spheroid self-gravity, 

Lsrf = LEddMsRF/MBH = Aa^c/{Gfg). (38) 

4.1. AGN triggering of star formation 

Negative feedba ck helps account for the black hole 
mass-<7 correlation (|Di Matteo et al. l l2008D and for the lu- 
minosity function o f massive galaxies ( Bower et al.ll2006l : 
ICroton et al.l [200^. More physics must be added how- 
ever to account for downsizing and efhcient star forma- 
tion in massive galaxies. The key may be AGN outflows 
that can trigger star formation by compressing dense 
clouds. These would precede the outflow phase which 
in this case is due to the combined effect of AGN outflow 
and triggered SNe. A prior phase of positive feedback is a 
possible new ingredient in feedback modeling and is mo- 
tivated by evidence ( admit tedly sparse but compelling: 
see e.g. iFeain et all (I200I ') for AGN trig gering of star 
formation. 

The following model is necessarily schematic pend- 
ing fully three-dimensional simulations of jet propaga- 
tion into a clumpy proto-galactic interstellar medium. 
We speculate that the triggering works as follows. Jet 
propagation into a clumpy medium de velops into an 
expanding, over-pressurized coc oon (Sax ton et 311120051 : 
lAntonuccio-Delogu fc Silkll2008f ) . This adds a potentially 
large multiplier to the efficacy of the BH-driven outflow 
for the following reason. The jet-driven plasma-flUed ra- 
dio lobe drives a cocoon that expands into the hot viri- 
alized component of the proto-galaxy at a speed Vco that 
is much larger than the velocity fleld associated with the 
gravitational potential well. Proto-galactic clouds that 
are above or near the Jeans mass will be induced to 
collapse. The cloud over-pressuring and resulting trig- 
gered star formation occurs at a rate much larger, by 1 
to 2 orders of magnitude, than is associated with nor- 
mal gravitationally-driven fueling, as would be appropri- 
ate to a star-forming disk galaxy. We show below that 
the star formation rate enhancement amount to a factor 

Vco/<^- Hence central massive black holes that have not 
yet grown by accretion to the limiting mass controlled by 
radiation outflow should still have a considerable impact 
on the evolution of the protogalaxy core. Incidentally, 
this early phase of black hole feedback makes observa- 
tional confirmation difficult, as discussed below. 

Numerical studies of radiative shock-induced cloud col- 
lapse reveal the cornplex interplay of h ot and cold gas 
(jMellema et al.ll2002l : lEYagile et al.ll2004l ). Here we focus 
on the implications for star formation via analytic con- 
siderations. For a simple ansatz, suppose that the trig- 
gered star formation occurs over the cocoon propagation 
time, tco that is much less than the dynamical time, td, 
tco ^ td ■ We show below that jet triggering of the co- 



coon expansion gives a star formation rate enhancement 
factor td/tco ~ Wco/fT- 

Although we shall generally assume Vco is a para- 
menter of the wind, consider the case in which the in- 
jected momentum flux is L/c for a wind that is mechani- 
cally driven but is originated via radiation pressure. We 
expect the mechanical jet luminosity to be a fraction 
Vw/c ~ 0.1 of the luminosity. More precisely, the ra- 
tio {Lmechlvw)/{L/c) is cqual to the optical depth r, 
due to a combination of Thomson scattering, line opacity 
and/or dust. To estimate r, an effective Rosseland mean 
opacity can be defined from these opacity contributions. 
This is valid if the jet is radiation pressure-driven. Now 



/ tAgn\ 1/2 
V Lsrf J 



(39) 



where we have set pg — /2TiGr^ . The previously de- 
rived star formation rate can be generalized for round 
systems with velocity dispersion a to 



{esNl(T)fcfciMg{GpgY'\ 



(40) 



Thus if the gas pressure is replaced by the AGN-driven 
pressure, we can see how a central AGN can boost the 
star formation rate by writing the pressure as 



PAGN = Lagn'^t^VcoT 



(41) 



The star formation rate is now given generally using 
equations (gOl) and (gl]) by 



AGN 



fcfcl^SN 



td 



r -1/2 

J a 



(— 

\Vagn 



1/2 



/ Lagn\ 



M. 



AGN 



fcfcl^SN 



M„ 



td 



1/4^-1/4 



1/4 / L 



-'AGN 



Lsrf 
(43) 

The AGN-driven enhancement factor is {pagn /PgY^'^ ~ 
(uco/o')t^/^. Note that esN c< a, so that the star for- 
mation efficiency coefficient (fraction of stars formed per 
dynamical time) is boosted considerably for spheroids 
relative to disks. The AGN luminosity explicitly drives 
star formation. It is the AGN-triggered star formation 
multiplier rather than the AGN itself that drives the 
feedback. The boost effect is generally important in the 
innermost spheroid, and globally important for super- 
Eddington AGN luminosities. 

Note that we must be careful that the AGN pressure 
does not blow away the ISM completely. For a SNe- 
regulated ISM, the pressure relates to Q and includes 
details of the SNe bubble evolution. But when the AGN 
jet provides the pressure it will over-pressure the clouds 
and turn them into stars with some efficiency. We need 
to relate the star formation from the triggered clouds to 
the AGN accretion rate and this would then affect the jet 
luminosity and pressure. Essentially, we raise the exter- 
nal pressure (here using the AGN) on the star-forming 
clouds. This allows us to regulate Q for the SN bubbles. 



1/2 



V Lsrf J 
(42) 

Here we have left Vagn an arbitrary variable to be 
applied to the appropriate situation such as jet, wind, 
outflow etc. In the radiation-driven case discussed above 
Vagn = Vco, the cocoon velocity, and using equations [39] 
and l42l we flnd 



1/4 
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If the porosity is maintained to be constant by the trig- 
gering of massive star formation, we have achieved the 
desired self-regulation between AGN activity and star 
formation. If Q is given by other means (and is of or- 
der unity), then star formation- regulated feedback fol- 
lows naturally. To show this point explicitly in our for- 
mulation of AGN driven star formation we return to the 
pressure-regulated model of the structure of the ISM, and 
instead of calculating the final size of SNell-driven rem- 
nant bubbles by setting the remnant's velocity, Va equal 
to CTg we set the external pressure PgOg^ equal to the 
mechanical pressure driven by the AGN, L/ Airr^Vco- As- 
suming, as above, that p = ag^/Gr"^ with an associated 
pressure Pg = Gg'^ jGr^ and replacing this with the AGN 
pressure, vagn given by vagn = Lagn I^T^VAGN-r^ we 
find that the required transformation to replace gas pres- 
sure with AGN pressure is 



1/4 



Lagn 



fid J \ VaGN 

It then follows directly that 

Qe-'^ = fci^SNjid (^~~^^ 
and for Q << 1 wc find 

Vagn 



L 



Q ^ fcl^SNJid 

For Q » 1 



1/4 



SRJid 



l^SRJid 

Lagn 



^SRJid 

Lagn 



g-in 



fcl^SNJid 



fVAGN \ ( L 



^SR.Jid 

Lagn 



1/4 



3/7 



3/7 



3/7' 



(44) 



(45) 



(46) 



(47) 



with a good fit to the range Q << 1 up to Q >> 1 given 

by 



g~in 



1 + fcl^SNJid 



/TA \ 1/4 / r X 3/7' 

/ VAGN \ I ^SRJid 



V 



Lagn 



(48) 



4.2. Cocoon Overpressure and the Bonnor-Ebert 
Condition 

Following the sim ple cocoon model of 
iBegelman fc"Cioffil (flOSl . we examine the effect of 
the power flow in the canonical two bidirectional jets 
diverted into a small nuclear cocoon and then via the 
cocoon pressure acting back on the central nucleus 
where gravitational instability and enhanced collapse 
and accretion may occur if the Bonnor-Ebert critical 
pressure is reached. For an isothermal gas distribution 
with velocity dispersion a as used throughout this paper, 
a jet opening angle 0j, an approximately ellipsoidal 
cocoon with axes a and h with a > b for simplicity 
powered by two thin jets with luminosity, Lj, we find 
that the cocoon pressure is given by: 

1/2 



P — 
^ r.n 



G 



(49) 



and the ratio of the cocoon pressure to the gas pressure 
in the central region (putting a ^ b) for simplicity is 

1/2 

(50) 



CO 



Pg \ LsRF J 



Now for the central isothermal core of gas pressured by 
the cocoon pressure, the ratio of the Bonnor-Ebert mass 
to the core mass is given by 



Mbe 
M 



1.18 



-1/4 1 // 



Lsrf) 



Thus, if vj ^ 0.1c, the core can be overpressured by 
the cocoon if 8j > 10(cr/c) ~ 10^^ when Lj ^ Lsrf- 
This may be another way to look at the feedback procees 
involving the growth of black holes along the Mbh — 
a line since at the Eddington luminosity Lj Lsrf 
implies M ~ Msrf- 

4.3. AGN winds 
The global mass loss for AGN-driven wind is given by 



M; 



\gal _ 



La 



GN 



out 



^ AGN 



(52) 



and for the radiation driven AGN-generated wind case 



Lag 



N 



G 



.-1/2 



(x 



GN 



SRF 



1/2 



(53) 



The outflow rate is proportional to the spheroid veloc- 
ity dispersion and to the square root of the AGN lu- 
minosity. The scaling of the outflow rate with regard 
to AGN luminosity implies that outflows saturate. It 
may be compared with observations of broadened fea- 
tures that demonstrate the presence of massive winds in 
ultraluminous star-forming infrared and radio galaxies. 

Making use of the AGN-enhanced star formation rate, 
we can express the outflow as a ratio using equations (|42p 
and ( [53|) generally as: 



MAGN 



n 1/2 

J 9 



f Lagn\^^^ f c 



fcfcl&SN{c/a) I \ LsRF J \Vagn 



3/2 



(54) 

It follows that the outflow rate from AGN is always of 
order the AGN-boosted star formation rate for the proto- 
spheroid. We may also compare the star formation- 
boosted global outflow rate with the AGN outflow. The 
AGN mass outflow rate is rjcMacc/ (vout), with vagn ~ 
0.1c. Hence it is of order the SMBH accretion rate. As 
expected, the mass flux associated with the galaxy out- 
flow dominates that from the AGN. The mass flux ratio 



IS 



Hence 



M, 



gal 



Mac 



c \ / Lsrf \ 



1/2 



(T ) \ Lagn ) 



(55) 



- 100 for AGN at the Eddington 
luminosity and n ^ IQ^gt /rUp. In order to allow for 
dust, if a factor r^^ is incorporated into the definition 
of Lsrf 7 this ratio is seen to be inversely proportional 
to the square root of the adopted (dust) opacity. The 
momentum flux ratio is 



MaccVAGN 



= Vfu 



c f Lsrf \ 
VAGN \ Lagn J 



1/2 



(56) 



We see that the momentum ejected from the AGN dom- 
inates over that in the global outfiow by a factor of a few 
for AGN near the Eddington luminosity Lagn Lsrf- 
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4.4. Downsizing 

The piston model enables downsizing of AGN and 
spheroids by coupling their growth. For some fiducial 
AGN energy conversion efficiency t] (~ 0.1), we note that 
j^AGN _ jjfP'M^^^ is a, measure of the BH accretion rate. 
Since L^^'w controls the star formation rate and is itself 
controlled by the black hole accretion rate, we infer that 
black hole growth and star formation triggering down- 
size together, provided Q is approximately constant due 
to AGN triggering of SN. The AGN driving of star for- 
mation overcomes the pressure suppression of porosity in 
the absence of the AGN. A large porosity also results in 
a wind. The required turbulent velocity field controls the 
accretion rate and might be specified by other physics, 
such as a merger, or even be due to the AGN itself. Let 
us try to make these assertions more quantitative. 

The AGN is the ultimate driver of the porosity. We 
need to connect AGN-induced star formation and out- 
flows to the black hole growth rate via the AGN lumi- 
nosity. The global outflow rate is 

= QhesNMg/td. (57) 
By momentum conservation, this must equal the global 
AGN-boosted outflow rate L^'^^/(ctc)- 

4.4.1. Downsizing for porosity-regulated star formation 

Incorporating the effects of porosity-driven star forma- 
tion means that the outflows must satisfy 

Mil = Q'fLia,/a^,,r-'f-'/'al/G. (58) 

The AGN luminosity is controlled by the accretion rate 
onto the central black hole. Mace- Our next step is to 
evaluate the black hole growth rate, Mace- This is the 
key to explaining downsizing. 

To repr oduce the downsizing phenomenon, observed 
for AGN fHasinge r et al.l 2005 ) and their massive host 
galaxies (Kriek et al .1 12007 ) to occur almost coevally, we 
need to understand why massive SMBH and spheroids 
form before their less massive counterparts. The re- 
quired scaling for L^'^^ or Maccr is reminiscent of the 
scaling found for proto-stellar jets. The magnetically- 
regulated disk phenomenon plausibly obeys a universal 
scaling law, that could equally apply to jets and out- 
flows from disks around SMBH. The proto-stellar scal- 
ingjs (|Mohantv et al.l l2005f ) Mace oc M^. lAllen et all 
()2006[ ) find that for the black holes that power the AGNs 
in massive ellipticals, the Bondi accretion rate is ap- 
proximately proportional to the jet power. The con- 
nection with outfiows and jets that are magnetically 
guided b y the wound-up field in the accretion disk pro- 
posed by iBanerjee fc Pudriti ()2006f ) is a generic scaling 
in their study of proto-stellar jets, M^^^ — fwMacc, 
with /u, ~ 0.1, for outflows associated with central ob- 
jects that range from brown dwarfs to super-massive 
black holes. 

The Bondi accretion formula therefore regulates 
SMBH growth and, implicitly, outflow. It yields Mbh = 
ttG^ {pg/a^) Mbh^ ■ A simple interpretation of this scal- 
ing is that for Bondi accretion, 

Monti U = Mace = ATT{GMla'' fpv (X (p/a3)Af^ (59) 
in combination with adiabatic compression, so that p cx 
. For the AGN case, we write Mace = ctM^j^ with 



a cx G^p/a^. Rewriting this we see that 
(IMbh j. -1 



dM^ 



Qhfu 



(60) 



Therefore if Q reaches a self-regulating constant and 
and /„, are also constant then the black hole and galaxy 
growth move together on a fixed trajectory in the Magor- 
rian plane as discussed below. It is this fixed trajectory 
that forces downsizing. 

4.4.2. Downsizing for SN energy injection 

Substituting further the relevant quantities for the case 
of SN energy input we find for the case without AGN 
feedback 



din Mbh _ . 1/2 fts\ _ 



where 



ts 



(61) 



(62) 



Kv^SNCfcfcl , 

Thus the critical parameter determining the logarithmic 
slope in the Magorrian plane is ^5. There is a critical den- 
sity Peril = {Gfgts)~^ above which black hole growth 
dominates and below which star formation dominates. 
This can be rewritten in terms of a critical velocity dis- 
persion if one takes Pcrit to be the density at the edge of 
the Bondi accretion sphere (the sphere of influence of the 
black hole, Rbh) namely p ^ (M/r^) ^ {MBH/rBH^} ^ 
{G~^a^MBH~^) giving then an equivalent 



M. 



SRF 



BH 



/■ 1/2 /. -1 /. -1 I Pel 
Jg fc fcl 



Pa 



1/2 



SN 



\msNVc 



(63) 

which is a satisfying combination of black hole-galaxy 
and ISM properties. Continuing to the case with AGN 
feedback, we flnd 



d In M BH 
din M* 



f 1/2 f ts\ ( VAGN \'^/'^ ( Ls 



RF 



Lagn 



1/2 



and in the radiation-driven case 
dlnMBH fi/i(ts\(<^\^'^ 



(64) 



dlnM, 



-1/4 



i- 



SRF 



ON 



1/4 



Corresponding expressions for the critical density and 
velocity dispersion in the AGN case can be readily ob- 
tained. At higher redshifts, galaxy systems can be 
denser, although much of the physics of the nuclear re- 
gions depends on local physics, and thus the dominant 
black hole growth phase may be more easily entered at 
higher redshift. The Bondi accretion formula can be 
rewritten as 

^ ^ MsRF \ 

Using the equation balancing inflow and outflow and 
the equation for the logarithmic slope in the Magorrian 
plane, we flnd 

/ M, \ (QfL\ ^ (t_s\ (VAGNy/^ ( Lsrf V'^ 

\MBHj\Pf^J \tdj\ c J [LagnJ ■ 

(67) 
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Eliminating (ts/td) we find 

f 9Il±^ ( ^bh \ _ fric\ / VAGN \ f Ledd \ f Msrf\ 
\I3UJ \ ts J ^\<j)\ c J [lagnJ V J 

and in this case we have used the momentum balance 
equation with momentum injection from SN balancing 
dissipation. For the radiation-driven case we find 

(QIlV (tmA ^ ^1/2 / Ledd^ ( Msrf V 

\PfwJ \ts J \LagnLsrfJ \ J ■ 

(69) 

In both the above cases there is no evidence for down- 
sizing even with constant Q. However, if there is an- 
other way to deduce Q for these turbulent AGN-driven 
multi-phase media and the momentum injection for the 
medium is taken up by the AGN then, as we shall show 
in the following, downsizing can occur naturally as a con- 
sequence of the turbulent ISM properties. 

4.5. Constraints on Evolutionary Tracks 

There is an obvious but powerful constraint on the be- 
havior of the evolutionary tracks in the observed mass 
\n[MBH] — ln[M] plane. We first emphasize this concept 
of tracks is implicit in our model and that there is a flow 
of points in the mass plane with evolutionary arrows all 
pointing in the direction of black hole growth. The slope 
p = d\n[MB h] / InldM] of the track cannot be negative 
since: (I) black holes can only grow in mass, and (II) 
galaxies only grow in mass (ignoring their small frac- 
tional mass loss). On dwarf galaxy scales the fractional 
mass loss can be easily incorporated but even there it is 
much less than of order unity. Therefore, for example, 
tracks cannot loop back from above the mean line with 
any slope less then zero after overshooting the mean line 
on a trajectory originating from below. Therefore any 
non-pathological track will spend most of its life on a 
track with a slope close to the mean. Thus we can as- 
sume p is approximately constant, Observationally p is 
of order unity. We now use this slope as a parameter in 
our time scale and evolutionary equations and find: 

fp\^f tBH \ _ /77£\ / VAGN \ f Ledd \ f MsRF \ ^ 

\(3J \tsj^\(j)\ c ) \LagnJ \Mbh J 
For the radiation-driven case we find 

fp\^f ^'Bh \ _ 1/2 f Ledd^ f Msrf \^ 

\P) \tsj~'' \LagnLsrf) \MbhJ 

(71) 

4.6. Why the AI^ dependence of accretion? 

We see that the parameter controlling accretion a, 
which is proportional to the phase space density, and also 
measures the specific entropy s of the initial gas distribu- 
tion, is specified by the physics of accretion. Specifically, 
a oc s~^/^, where s = kTn~^^^. In terms of a polytropic 
equation of state, a is constant for 7 = 5/3. 

4.6.1. Hot Phase Entropy 

In any dissipative multi-phase medium, the entropy 
cannot be strictly constant. However to demonstrate 
that entropy is indeed slowly varying in the hot phase 
of the system system, consider thermal balance between 



gravitational accretion heating, which is also propor- 
2 tional to the resulting power in the AGN outfiow and 
/rjfAssociated heating, and atomic cooling. One obtains 
civ^/G — p'^r^\{v), where \{v) represents the cooling 
rate per atom per unit density. In the range of interest 
where hot gas dominates the gas pressure on galaxy and 
cluster scales, X{v) is weakly varying, e.g. X{v) oc v for 
thermal bremsstrahlung cooling at T ^ lO^K and is ap- 
proximately constant over 10^ - lO'^K. In fact to avoid 
fragmentation, a n ecessary condition for effec tive central 
black hole growth (|Lodato fc Nataraianll2006[ ). one needs 
to be at r > lO^K and have 7 > 4/3. Above lO'^K, ap- 
propriate to massive galaxies and clusters, s cx T^/^ is 
found to be slowly varying and this helps account for 
the central entropy "floor" in clusters. This results in 
the maximum accretion rate being in the core. Black 
hole growth by Bondi accretion indeed requires constant 
specific entropy flow, which explains why Mace oc Mg^. 
These arguments should apply on massive galaxy scales 
where the gas pressure is controlled by the ISM hot gas 
phase at the outer boundary of the flow. 

4.6.2. Cold Phase Entropy 

Similar arguments apply to a cold phase in a multi- 
phase medium. Here the cold clouds themselves, envis- 
aged as bound self-gravitating entities that move on bal- 
listic orbits, act like massive particles whose dynamics 
can be described by Bondi accretion. For supersonic tur- 
bulence there are three points to note: 

(1) the system is highly dissipative, and so momentum 
conservation not energy conservation is the rule; 

(2) the density as described in the PDF is essentially di- 
mensionless and only measured in units of the square of 
the Mach number (M)^; 

(3) continuous energy and momentum input comes for 
the central source so that the PDF structure of the 
medium remains statistically robust. This means there 
will be the same number of clouds with the same mass 
function even though there is continuous creation and 
destruction. If a cloud is dissipated then another cloud 
will be created to take its place in the ensemble. 

Finally, the dependence discussed above occurs 
\ ' ^iover a very wide range of sources and environments both 
relativistic and non-relativistic and with power sources 
ranging from proto-stars to micro-quasars to quasars. A 
clear invariant seems to be supersonic turbulence gener- 
ated by a central source. For the cold supersonic phase 
p ^ and a ^ M and thus we expect the phase space 
density be p/cr^ ~ A^^^. Note that over a wide range 
of the systems discussed from proto-star to quasars the 
value of 3 — 10 is appropriate. In addition the cool 
gas over a wide range of conditions is at similar tempera- 
tures. The cool gas is carried in packets (in clouds) whose 
number density and mass function are properties of the 
supersonic turbulence and the Mach number. Thus over 
this wide dynamic range the variation in the parameter 
a may not be large. 

4.6.3. Multiple phases 

For cloud populations with different properties their 
effective temperature is associated with the velocity dis- 
persion of the cloud ensemble, as opposed to the gas 

kinetic temperature. In this case, s = "^Miafn^ 
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where we sura over hot and cold cloud components, and 
the Bondi accretion rate for a two-phase medium is now 



4.7. From Bondi accretion to star formation 

We now develop the interplay between the Bondi ac- 
cretion rate parameter a and the porosity-driven star 
formation rate. We show that constant Q implies con- 
stant a, and vice versa. This is the key to under- 
standing coupled downsizing for spheroids and super- 
massive black holes. Define the black hole growth time 
by tBH = MBH/Macc = l/(aMsff). 

4.7.1. Significant downsizing for porosity-regulated star 
formation 

Using the above equations ([TO]) and ([55]) which relate 
Mout and Macc^ we find that 



Q'fL 

fw 

f Mbh \ 
[msrfJ 



which becomes 

fw 

Mbh 



yc \ / tBH 

CTfldJ \ ts 



M, 



BH 



Mn 



31/7 



(72) 



rjc 



ft_BH\ 
CTftdJ \ ts J 



MsRF 



/ 1\ T \ / T \ 31/28 



V Ls 



RF 



(73) 



Alternatively using the p-parameterization discussed in 
section l475l we find 



{QpY 

M„ 



^ fid J \ ts J 

( Mbh \ /f7^d\3i/7 



, . , > , • (74) 

^MsRF J \MsRF, 

Constant porosity therefore guarantees downsizing, 
since 

m 

For Eddington- limited accretion, Mbh — -Wsi?F -The 
ratio of black hole growth time to Salpeter time decreases 
with increasing black hole mass at constant Q. Notice 
also that 



^-31/28 31/28 
°^ ^SRF °^ ^^^BH 



M 



BH- 



(75) 



l/a = tBH Mbh oc Q 



-4»,r2 ,^-59/28 ,-—3/28 



^BH^SRF ^"^BH 

Hence constant porosity also favors Bondi accretion since 
a ~ constant. Since a ~ lOaftd for a massive spheroid, 
we also infer that for Eddington-limited accretion, the 
porosity must be of order 10 percent if the black hole 
growth time is of order the Salpeter time. Indeed, we can 
equate these time-scales and infer that porosity depends 

weakly on black hole mass, Q ~ M^^f . In the absence 
of AGN feedback, the derived star formation law yields a 
ratio of star-formation time-scale to dynamical time that 
is proportional to ag and hence approximately constant, 
independently of galactic mass. This suggests that in 
starbursts when AGN play no role, there should be no 
downsizing, as would be expected if internal processes 
such as those associated with formation of massive star 
clusters were to dominate. 



Also, the Eddington ratio can be written as 
fEdd = — = vc ( I Mbh- 



UtOL 

rripG 



LEdd 

Hence constant a is consistent with the obse rved trend 
meas ured in the Eddington ratio ( Alonso-Hcr rero et al.l 
120081 ). The Eddington ratio is found to be lower for AGN 
than for QSOs due to a combination of reduced host 
galaxy (and SMBH) mass as well as AGN feeding. 

4.8. Why does porosity self-regulate? 

First, we give a qualitative argument for self- 
regulation. If the porosity Q is low, the jet is blocked, and 
the turbulent pressure is enhanced. Weak shocks prop- 
agate ahead of the cocoon and squeeze self-gravitating 
clouds over the Bonnor-Ebcrt stability limit. Star for- 
mation is triggered and the resulting supernovae drive up 
the pressure. Blow-out most likely occurs of the residual 
gas. This in turn increases Q ~ 1. However dense clouds 
can now fall in unimpeded by intervening dense cold gas. 
The infall reduces Q, drives accretion and resurrects a 
strong jet. 

Now, it is reasonable to assume that the cold phase 
is defined by a minimum density ricr set by dissipative 
cooling and molecular formation, and that is not strongly 
dependent on met allicity and ionization fra ction either 
in the molecular ( Norman fc Spaani Il997f ) or atomic 
(jWolfire et al.l l2003l ) phases! The shape of the density 
PDF is well represente d by a log normal distribution 
(jWada fc Normanll2007l ) with two parameters, dispersion 
of the distribution and its amplitude. When calculating 
the filling factor of either hot or cold gas, one integrates 
the lognormal PDF below or above ricr- At fixed n^n the 
filling factor of either hot or cold gas depends only on one 
parameter, the dispersion. In general, for the isothermal 
case considered here, the dispersio n is a linear function 
of In Al, for high Mach numbers (jKrum holz fc McKed 
l2005f ). In the adiabatic case the lognormal PDF disper- 
sion is independent of the Mach number. The lognormal 
form is retained for supe rsonic turbulence both w ith and 
without star formation (jWada fc Normaiil l2007f ) . Even 
with significant feedback and increased Mach number, 
Q only varies logarithmically with Mach number, and 
therefore any change of filling factor in systems with de- 
veloped supersonic turbulence is relatively slow as the 
turbulence is increased. In summary, the volume filling 
fractions of cold gas e~'^ and hot phase 1 — e~'^ depend 
only logarithmically on the Mach number. Hence the 
porosity is plausibly constant. Quantifying this argu- 
ment we define the hot-phase filling factor in a turbulent 
supersonic medium to be the volume of the turbulent 
medium with a lower density than the mean by a factor 
Vh < 1. Using the usual lognormal probability density 
function (PDF) 

/(n) = ^=^exph^i^] (76) 

where for supersonic turbulence we use 

at^ = In [1 + XtM^] (77) 
with M being the Mach number, and fo r the parame- 
ter At we use At = 3/4 f Krumholz fc McKee. ,20051 for 
numerical estimates. We define 

^ no [1 + X,M^] ^''^ 
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where hq is a reference density, and we find an expression 
for fh 



fh = erfc 



In 


j_ 







V2, 



CTt 



(79) 



where often the first term of the asymptotic expansion 
for erfc[x] = exp [—x'^J/x^/tt + .... is a useful guide. For 
the Mach number, Ai = 3 — 10 and for the under density 
parameter of the hot phase relative to the mean we use 
Vfi = 0.1 and find the hot phase filling factor is ~ 10 — 
20%. As the Mach number increases, fh increases in turn 
as the width of the PDF increases. Since fh = 1 — exp ~Q 
we find that 

Q = In [1 - fh] (80) 
which for Q << 1 becomes 

Q = In [1 + fh] (81) 



givmg 



Q-ln 



1 + erfc 



In 



V2< 



(82) 



This shows how Q increases as the Mach number in- 
creases. Proceeding further, we can generalize our pre- 
vious expression for Q, including both the competing ef- 
fects of the star formation and supersonic turbulence, to: 



Q-ln 







■ 1 ■ 




1 -f erfc 


In 


V 




.V2CTt. 





+ fclf-SN.Jid 



<yfid 



12/7' 



(83) 



Feedback from AGN is now readily incorporated by sub- 
stituting Ai — {p jo jid)^ and then making the substitute, 
as before, for the multi-phase medium under the action 
of the mechanical and radiative pressure originating in 
the central source giving 



1/2 



Lagn 



1/2 



\VaGN J \LsRJid 

and thus using equations ([77)1 and ([M]) we obtain 



(84) 



<^t,AGN 



In 



l + Xt 



^AGN 



\Vagn J \LsB., 



fid 



(85) 

In the AGN case it is now clearly quantified how the AGN 
pressure is reducing fh by confining the hot SNR bub- 
bles but on the other hand how the AGN increases the 
Mach number of the turbulence and therefore broadens 
the PDF distribution, thus increasing fh- 

4.9. Star formation rate 

Global gas consumption is dominated by star forma- 
tion, and locally by SMBH growth. The two are con- 
nected via the piston model and suggest a possible self- 
regulation loop for both spheroid and SMBH growth. We 
now show that the time sequence underlying the Magor- 
rian relation can be interpreted in terms of the ratio of 
spheroid to SMBH growth rates. 

The star formation (or gas consumption) rate is from 
equation (j42|l 



1 _ ^SN f Lagnt 

t* td V LsRF 



1/2 



(86) 



/ Lagn\ 

\ LsRF J 



1/2 



1/2 



(87) 
(88) 



where esN = £sn/<^ is a constant depending only on su- 
pernova properties and the IMF. This controls spheroid 
growth and demonstrates spheroid downsizing via the 
AGN-enhanced pressure and associated star formation 
rate. 

Writing the AGN luminosity as Lagn = V^bhc^, the 
Eddington luminosity can be expressed in terms of the 
Salpeter time as tjMbhc^ /ts- We now rewrite the star 
formation rate expression and obtain 



BH 



t* J XesNcrTVc^t 



(89) 



The predicted normalization of the Magorrian relation 
agrees with the local value and slope for canonical pa- 
rameter values (Stot ~ S*, 77 ~ 0.001, esN ~ 0.1). Of 
course, there is considerable uncertainty due to possible 
variations in the initial mass function, supernova energy 
and star formation timescale. 

In fact, the relevant SMBH measure in distant objects 
is Lagn rather than Lsdd- Let us make use of Lsrf 
as a fiducial luminosity, in effect a proxy for cr^.We now 
rewrite the preceding expressions to obtain 

Lagn / Lsrf = ^sNifg^d/t*)"^ ~ {td/t^)"^. (90) 

We also have for the AGN-boosted star formation rate 

M^"^"" ~ esNMgQ {LAGN/LsRFf'^ Mg^. (91) 

This is of course the optimal rate. We also see that 
j^j-AGN oc This is not inconsistent with the observed 
dependence of outflow velocity on star formation rate 
(Martin 2005). 

We may consider the case of a recently detected kilo- 
parsec scale starburst at z = 6.24, hosted by a lu- 
minous quasar which has spatially resolved [CII] emis- 
sion as well a s a large resery oir of CO-detected molec- 
ular gas ( Walter et al. Il2"009f ). Other similar high z ob- 
j ects, detected in CO, are believed to be super-Magorrian 
(jMaiolino et al.ll2007t l. This quasar host galaxy also has 
a star formation rate of ^ 1000MQyear~ikpc~2, an order 
of magnitude larger than is typical of starbursts without 
luminous AGN. For comparison, Arp 220, a low redshift 
starburst hosting an AGN of luminosity comparable to 
its starburst, has a similar surface brightness in star for- 
mation but only over a 100 pc scale. It is tempting to 
infer, admittedly with only two well-mapped examples, 
that we may be viewing AGN boosting of star forma- 
tion, with the phenomenon being greatly magnified at 
high redshift for the most massive objects. 

We infer that if the SMBH mass is super-Magorrian, 
then < ts and t^, < td, and star formation is very 
efficient. This seems to be the case at high redshift. 
The preponderance of data indeed suggests that the lo- 
cal relation become s super-Magorrian prior to z ~ 2 
(|McLure et al.l 120061: IWoo et al 1 120081) and indeed per- 
sists to z ^ 5 (jMaiolino et al.l 120071 ). Coevolution of 
AGN accretion and the co-moving star formation rate 
densities occurs to z ^ 2, but the accretion rate fa lls off 
relatively towards higher redshift (jSilverman et al.l2008 ) . 
Comparison of the cosmic star formation history and 
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AGN accretion rates in comoving number density as a 
function of luminosity suggests that the peak in mas- 
sive black hole growth rate occurs several Gyr prior to 
the star formation peak and that d ownsizing at z < 1 
is due to diminishing accretion rates (jBabic et al.ll2007l ). 
Sub-millimeter galaxies (SMGs) are an exception. SMGs 
at z ~ 2 contain SMBH that are under-massive rela - 
tive to the Magorrian relation (|Alexander et al.ll2008D . 
This is suggestive of triggered star formation, which re- 
duces f*, and may be appropriate in major mergers that 
generate dense central gas environments where porosity 
feedback is suppressed. Note also that the peak in the 
major merger r ate also precedes the peak in cosmic star 
formation rate (jRvan et al.ll2008l ). and is approximately 
consistent with the peak in comoving AGN accretion rate 
density. 



5. SUMMARY AND CONCLUSIONS 

In summary, a general and robust treatment of disk 
star formation is developed from a cloud collision model. 
The Schmidt-Kcnnicutt law emerges naturally for star- 
forming disks. We predict that there is an inverse re- 
lation between TuUy-Fisher law and Schmidt-Kennicutt 
law residuals. A multi-phase treatment of supernova 
feedback leads to a turbulent pressure-regulated general- 
ization of the star formation law that is applicable to 
gas- rich starbursts. Negative feedback from star for- 
mation occurs in disks under turbulent pressure regu- 
lation. In combination with a cloud collision model, the 
Schmidt-Kennicutt law can be understood in diverse en- 
vironments, spanning quiescent disks and starbursts. En- 
hanced pressure, as expected in merger-induced star for- 
mation, enhances star formation efhciency. An upper 
limit is derived for the disk star formation rate in star- 
bursts that depends only on the IMF and on the ratio 
of global to cloud pressures. For clouds in approximate 
pressure with interstellar medium and a local IMF, we 
infer a limiting gas surface density of ~ IOOOMq pc~^. 

We extend these considerations to the case where the 
interstellar gas pressure in the inner galaxy is dominated 
by outflows from a central AGN. The star formation rate 
is pressure-driven and depends on the excess pressure 
applied by the AGN outflows. During massive spheroid 
formation, AGN-driven winds trigger star formation, re- 
sulting in enhanced supernova feedback and outflows. 
Downsizing is predicted to be a consequence of AGN- 
driven positive feedback. Our most important results 
refer to downsizing, for which we provide a new interpre- 
tation in terms of Bondi accretion feeding of the central 
black hole. 

The specific accretion rate is proportional to the black 
hole mass. We found that Bondi accretion results in 
Mbh oc CT^^^'^Q^{tBH/ts)- This means that if poros- 
ity self-regulates to be constant, black hole growth pro- 
ceeds rapidly until it saturates at the Magorrian relation 
Mbh (t'* due to blow-out. Black hole downsizing oc- 
curs if a is approximately constant. We clarify this as 
follows. 

There are three specific rates that define our model. 
The Salpeter rate tg^ is constant, the black hole growth 
rateisl/<BH = (l/t5)(-LyiGA'/-^£;rf(i), and the star forma- 



{GLAGNTy^^ca^fgr'^'. Hence 



1/2 



^t1/2 (92) 
td 



tion rate is 1/t* = esN 

tsH _ f-SNts 
t* fgtd \LaGnLsrf 

and we infer that 

Lagn/Lsrf 

This shows that the black hole growth rate and star for- 
mation rate are coupled. At given a, there is a critical 
AGN luminosity, above which AGN-triggered star for- 
mation rates dominates over the black hole growth rate. 
This critical luminosity is found to be 



[td/Uf ~ (ts/tBH? 



(93) 



L 



Edd 



fg td M a \ J V / 



(94) 



At super-Eddington luminosities, AGN-triggered star 
formation dominates. 

The model contains two characteristic luminosities 
which are functions of a. The Eddington luminosity, 
if combined with the quenching assumption, scales as 
a*. The AGN-triggered star formation luminosity is 
eic^^SNcr'^{LAGNT)^^'^{Gcfg)~^/'^. This is proportional to 
a^/^Msijcr^. Adopting Mbh oc ct* and a cx cr'*/'^ if Q is 
constant, we find that L^'^^ cx a^^^"^. Hence we again 
infer a critical Eddington luminosity above which trig- 
gered star formation dominates the luminosity of the 
system. This guarantees efficient star formation for lu- 
minous AGN. Moreover if Mbh/o''^ increases with in- 
creasing redshift, as inferred from the downsizing of the 
black hole growth rate, spheroid star formation down- 
sizes both in mass and in efficiency. 

The ratio of AGN to star formation luminosity is 



■qM^ 



bh 



Lagn 



dM, 



bh 



dM* 



(95) 



This yields an arrow of time in the form of tracks in the 
Magorrian diagram. If for each data point in the Magor- 
rian plane, defined by black hole and bulge mass, one 
can separate star formation and AGN luminosity, then 
the ratio gives a vector and hence an arrow of time. This 
is of course the instantaneous trajectory of points in the 
Magorrian plane as viewed at any given epoch. However 
statistically the vectors should provide the flow of galaxy 
bulges towards the Magorrian relation. It would be in- 
teresting to construct Magorrian flow diagrams binned 
over several redshift ranges. This would provide insight 
into the cosmological evolution of the flow of points in 
the Magorrian plane. 

AGN-enhanced star formation is mediated by turbu- 
lent pressure and relates spheroid star formation rate to 
black hole accretion rate. As the AGN pressure is in- 
creased, via jet/cocoon pressure-driven interactions with 
the ISM, the induced star formation rate is correspond- 
ingly boosted. Downsizing for spheroid formation is a 
consequence of massive black hole downsizing. The ob- 
served relation between black hole mass and spheroid ve- 
locity dispersion is obtained, with a coefficient (the ratio 
of Salpeter time to gas consumption time) that provides 
an arrow of time. Highly efficient, AGN-boosted star 
formation is favoured at high redshift. Outflows are of 
order the AGN-boosted star formation rate and saturate 
in the super-Eddington limit. 
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Wc end with some relevant questions for observers that 
are pertinent to our model. Is AGN activity correlated 
with the star formation rate? Was spheroid formation 
and black hole growth coeval and symbiotic? Which 
came first, if either? Is the reach of the AGN too lo- 
calized to globally affect star formation? Is or was feed- 
back significant in radio-quiet AGN? Are the youngest 
radio sources, notably the GPC sources, templates for 
the earliest stages of AGN feedback, and if so, is there 
associated triggering of star formation? If star formation 
is triggered by radio cocoons, where is the evidence for 
a post-starburst stellar population in old radio lobes? if 
the efficient mode of star formation is due to coherent 
cocoon triggering as argued above, what is the evidence 
for spatially and temporally coordinated episodes of star 
formation in well-studied examples such as the Anten- 
nae? Has the trigger of positive feedback in ultralumi- 
nous starbursts disappeared, due to a short duty cycle, 
or is it well and truly buried in embedded AGN nuclei? 
Are the associated outflow rates from AGN/starbursts of 
any significance for balancing the baryonic budget of the 
miiverse, and if so, where does the enriched debris end 
up? If the gas remains in the halos of massive galaxies, 
as essentially all current semi-analytic simulations pre- 
dict, why hasn't it been seen? And for the modellers 
(who are effectively observers of the computer), how can 
AGN feedback simulations possibly play any predictive 
role if one has to begin with massive seed black holes of 
dubious heritage and uncertain fate? Perhaps our ana- 
lytic discussion will motivate more realistic recipes and 
sub-grid physics prescriptions for future generations of 
feedback simulations. 

Many of these questions were inspired by discussions 
with participants at the Oxford-Catania Workshop in 
Vulcano, May 2008, on the Interface between Galaxy 
Formation and AGN. We gratefully acknowledge all of 
them as well as the role of Vincenzo Antonuccio-Delogu, 
who organized a brilliantly topical meeting in the unfor- 
gettable setting of the Aeolian Islands as well as provided 
an opportunity for us to complete this paper. 

The research of CN was funded in part by NASA grant 
GO 6-1730-X and a JHU/APL collaborative grant. 

6. APPENDIX:LIST OF SYMBOLS 
a Bondi accretion rate parameter: phase space density 

f3 factor for black hole growth timescale 

CsK Schmidt-Kennicutt law coefficient 

Ctf TuUy-Fisher law coefficient 

EsN kinetic energy of a SNell 

esN star formation eflaciency factor 

^SNjid star formation efliciency at fiducial velocity 

e; energy release per unit rest mass 

Cm mechanical energy release per unit rest mass 

e„ nuclear burning efliciency per unit rest mass 

esN modified star formation factor 

Eg gas mass surface density 



cTg gas velocity dispersion 

Ec/ cloud mass surface density 

fci gas fraction in clouds by mass 

fc cloud volume filling factor (= e^^) 

fL hot gas loading factor 

/; luminosity per unit mass of massive star formation 
f-u, fraction accreting material that flows out in wind 
H disk scale height of gas 

I-band luminosity 
Aeff effective cooling function for SN bubbles 
A cooling rate per atom per unit density 
LsRF self-regulating feedback luminosity 
Ledd Eddington lumimosity 

LsR.fid self-regulating feedback luminosity at flducial 

velocity dispersion 

Lmech mechanical luminosity 

Lj jet luminosity 

At numerical constant 

msN mass formed in stars per SNell 

Mout mass outflow rate 

MsRF self-regulating feedback mass 

Mbh black hole mass 

Mbe Bonnor-Ebert mass 

^out g^l^^y mass outflow rate 

]\/[AGN AGN-enhanced star formation rate 

M. turbulent Mach number 

Uh fiot phase underdensity parameter 

f2 angular rotation of the disk 

p logarithmic slope in the black-hole mass vs galaxy 
plane 

Pg gas pressure 

Pel pressure in a cloud 

Pco cocoon pressure 

Q porosity of the interstellar medium 

R radius of disk 

Ra maximum SN bubble radius 

Ro,to,vo, Co defined constants for SN bubbles 

Pg gas mass density 
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Sci covering factor of clouds 

E» mass surface density in stars 

Titot total mass surface density of the disk 

a fid fiducial velocity dispersion for SN driven ISM 

a velocity dispersion of system 

s specific entropy 

at PDF dispersion 

gas consumption time 
is modified Salpeter time 
tcoii cloud collision timescale 
ta time for SN bubble to reach maximum radius 
tc cooling timescale of SN bubble 
td dynamical timescale 
tco cocoon propagation timescale 



r Rosseland mean optical depth 

Ibh Black hole growth time 

ts Salpeter time scale 

0j jet opening angle 

Vc supernova velocity at strong cooling 

Vr rotation velocity of the disk 

Vco cocoon velocity 

Vyj wind velocity 

Vagn velocity of AGN outflow 

vj jet velocity 

X interstellar cloud pressure enhancement due to self- 
gravity 
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